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Abstract 

A phenomenological theory of luminescence properties of one-dimensional resonant photonic 
crystals is developed within the framework of classical Maxwell equations with fluctuating polar- 
ization terms representing non-coherent sources of emission. The theory is based on an effective 
general approach to determining linear response of these structures and takes into account forma- 
tion of polariton modes due to coherent radiative coupling between their constituting elements. 
The general results are applied to Bragg multiple-quantum-well structures, and theoretical lumi- 
nescence spectra of these systems are compared with experimental results. It is shown that the 
emission of such systems can be significantly influenced by deliberately introducing defect elements 
in the structure. The relation between absorption and luminescence spectra is also discussed. 
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I. INTRODUCTION 



A possibility to influence the emission of light by tailoring the dielectric environment 
of emitting objects has been attracting a great deal of attention since a pioneering work 
by Yablonovitch.- In Ref. Q it was suggested that a three-dimensional periodic modulation 
of the dielectric constant can result in formation of a photonic band structure, consisting 
of allowed and forbidden photonic bands in analogy with electronic band structure. One 
of the important consequences of the band structure is a modiflcation of electromagnetic 
density of states, which can be used to suppress or enhance the rate of the spontaneous 
emission of emitters embedded in a photonic crystal. While such drastic effects as the full 
inhibition of the spontaneous emission proved to be difficult to achieve,^^'- there is still a 
growing interest in emission properties of photonic crystals,-'^'^'^'^ which even in the absence 
of complete photonic band-gap may result in a significant modification of properties of 
emitted radiation. If one is not looking to achieve the full inhibition of the spontaneous 
emission, the systems in which a periodic modulation takes place in only two or even one 
dimension are also of interest, because even though they cannot confine light completely, 
they do modify emission patterns for particular directions, which can be useful for various 
applications. 

One-dimensional structures attract a particularly great attention, firstly, because they are 
easiest to manufacture, and secondly, because they allow for a detailed theoretical descrip- 
tion. These two circumstances make one-dimensional structures most suitable candidates for 
a number of applications that do not require modification of the electromagnetic properties 
in all three dimensions. At the same time, certain properties of one-dimensional structures 
are typical for two- and three-dimensional systems as well, and, therefore, these structures 
provide a convenient testing area for understanding some general properties of media with 
spatial modulation of the dielectric function. 

Most of the previous works addressed the problem of the spontaneous emission in ID 
photonic structures from the perspective of individual emitters embedded in special dielectric 
environments such as superlattices or a Fabry-Perot cavities (see, for instance, Refs. yQiiii, 
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and references therein). Theoretical analysis of these situations is based on the assumption 
that the structure of photonic modes is determined solely by the periodical modulation of the 
dielectric function. The interaction between the photonic modes and emitters is considered 
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in this approach as weak in a sense, that it does not affect the structure of the photonic 
modes, and can be treated perturbatively within the framework of Fermi Golden Rule. 

These assumptions, however, break down in the important case of so called reso- 
nant photonic crystals, that have begun attracting considerable attention in recent years. 
These structures are composed of periodically distributed structural elements containing 
dipole active internal excitations, which co-exist with periodic modulation of the dielectric 
constant.— 'i^'i^'i^'i^'i^'i^'22i2ii2^ Multiple-quantum- wells^^ (MQW) present one of the popu- 
lar realizations of such structures with the one-dimensional periodicity. Excitons confined 
within quantum wells provide optically active excitations, and the contrast between re- 
fractive indices of wells and barriers is responsible for periodic modulation of the dielectric 
constant. When the period of the structure satisfies a special, so called, Bragg condition, the 
interaction between light and excitons cannot be considered as weak, and cannot be treated 
with the help of the Fermi golden rule. Excitons and periodic modulation of the refractive 
index play in such structures equally important roles in the formation of the photon modes, 
which should be more appropriately called polaritons. The Bragg condition, in its most 
general form can be written as A0p(co'o) = vr, where A0p is the change of the phase of the 
propagating electromagnetic wave over one period of the structure calculated at the exciton 
frequency uj^r^. If one neglects the refractive index contrast (optical lattice approximation) 
this condition can be rewritten as cu^/d = tt,— where d is the period of the structure, and 
c is the speed of light in the medium. 

The emission of light in such structures differs significantly from the cases considered 
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in Refs. I9l jl0ljllljl2l . The emitters of light in resonant photonic crystals affect the spatial 
structure of electromagnetic modes as much as the modulation of the refractive index. As a 
result, the processes of light emission by a particular quantum well and its propagation inside 
the structure should be considered on equal footing. To develop a theoretical formalism for 
dealing with such situations is the main objective of this paper. While focusing on the 
luminescent properties of Bragg MQW structures, we will treat them in a broad context 
of resonant photonic crystals. This allows us to develop a universal theoretical formalism, 
applicable to essentially any type of one dimensional structures with periodically distributed 
emitters. 

Since we will be interested in the effects of photonic environment on the luminescence 
rather than in a microscopical description of the processes of the exciton relaxation and 
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recombination, we will base our theory on macroscopical Maxwell equations with a non- 
coherent polarization source term, which would simulate non-coherent exciton population 
created by a non-resonant pumping. The phenomenological nature of our work distinguishes 



it from recent paper 



25l . in which microscopic theory of spontaneous emission of Bragg 



MQW structures was based on an approach, in which both electron-hole dynamics and 
electromagnetic field were treated quantum-mechanically^. Phenomenological nature of our 
theory allows for establishing a direct and clear connection between global optical properties 
of MQW structures and their luminescence spectra, which is important for understanding 
an intimate relationship between geometrical structure of MQW systems and their emission 
properties. On the other hand, properties of the polarization source entering our calculations 
as a phenomenological object can only be established either from independent experiments 
or on the basis of microscopical theories of the kind developed in Ref. |25. 

The theoretical formalism developed in this paper solves also a more general problem of 
calculating an optical linear response of finite ID resonant photonic crystals. Usually, the 
linear response is studied using Green's function formalism, based on spectral representation 
of the Green's function. The later, however, is not very well suited to deal with finite struc- 
tures whose normal modes cannot be considered as eigenfunctions of a Hermitian operator. 
In our approach we develop a method of relating Green's function of the finite resonant pho- 
tonic crystal to transfer matrices describing reflection and transmission properties of these 
structures. 

The paper has the following structure. In Section [III we formulate the basic equations 
describing the distribution of the electric field in the structure and introduce the basic ele- 
ments of our transfer matrix approach. In Section IIIII we conclude the presentation of the 
general formalism by deriving a general expression for the respective Green's function. In 
Section [IV] the general formalism is applied to the problem of the exciton luminescence spec- 
trum in resonant photonic crystals. In addition to considering an ideal periodic structure, 
we also discuss the role of inhomogeneous broadening, modifications in the luminescence 
induced by a deliberate introduction defects in the otherwise periodic structure, and rela- 
tion between luminescence and absorption spectra. The paper is concluded with Appendix, 
where we discuss relation between our phenomenological approach and microscopic theories. 



4 



II. GENERAL SOLUTION OF MAXWELL EQUATIONS IN ONE- 
DIMENSIONAL RESONANT PHOTONIC CRYSTAL WITH POLARIZATION 
SOURCES 



A. Maxwell equations for multiple-quantum-well system with sources of polariza- 
tion 

Our approach to description of emission properties of MQW structures is based on solu- 
tion of classical Maxwell equations for monochromatic field with frequency uj of the following 
form: 

V X V X E = — [n^{z)'E + 47rPexc + 47rF] , (1) 

where the coordinate z is chosen to represent the growth direction of the structure, and 
n{z) is the periodically modulated background index of refraction: n{z + d) = n{z). The 
polarization in this equation is presented as a sum of two terms. First of them, Pexc, 
originates from optical transitions, with frequencies within the spectral region of interest and 
has, therefore resonant behavior. The second contribution, F(r), arises due to other emitting 
transitions, which are non-resonant at given frequencies. In typical experimental situation 
involving /// — V quantum wells, the resonant term corresponds to polarization due to Is 
heavy-hole excitons, and the non-resonant contribution can be considered as a contribution 
from all other optically active transitions between states of electron-hole system. Such a 
separation of polarization into contributions from different optical transitions is possible if 
one neglects Coulomb correlations between exciton states and electron-hole plasma. These 
correlations are important for for highly excited states of semiconductors^^, but can be 
neglected if concentration of photo-excited electron-hole pairs is not too high, which is the 
situation considered in this work. 

Semi-classical description of light-exciton interaction in a single quantum well (see, for 



instance, Refj27l) indicates that the exciton contribution to polarization of a single QW can 
be presented in the following form: 



,p) + Sr„(p) 



(2) 



Here index m numerates quantum wells, p is two-dimensional position vector perpendicular 
to the growth direction of the structure, and the wave function of the exciton localized in 
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the m-th well, is taken in the form ^m{z) = $(-2 — Zm), where Zm is the position 

of the center of the m-th well. The expression proportional to the electric field in this 
equation describes direct optical excitation of excitons, while the second term, represented by 
function Sm(p), introduces an additional source of exciton polarization due to non-radiative 
processes. Depending on the properties of this term it can describe either coherent or non- 
coherent emission. More detailed description of function 'Sm{p) in relation to the description 
of luminescence, as well as general justification of the phenomenological approach to this 
problem is given below in Section HVl Here we will treat Sm(p) as an arbitrary function 
responsible for generation of exciton polarization. 

Writing Eq. ([2]) we explicitly take into account that the dipole moment of heavy-hole 
excitons is oriented in the plane of the well and, therefore, only components of the field 
perpendicular to the growth direction, E^, enter the expression for coherent exciton po- 
larization. The intensity of the exciton-light interaction is characterized by the exciton 
susceptibility Xm{^)- Neglecting the exciton dispersion in the plane of the quantum well 
and the inhomogeneous broadening the susceptibility can be written in the form 



where a is the exciton-light coupling parameter proportional to the exciton dipole moment, 
ujmo is the exciton resonance frequency in the m-th well, and 7 is the homogeneous broad- 
ening of the exciton line. 

The exciton polarization of the entire MQW system is the sum of polarizations of indi- 
vidual wells: 



Here we assume that the period of the spatial arrangement of the quantum wells coincides 
with the period of the modulation of the dielectric function, such that the distance between 
adjacent wells is z^+i — z^ = d. In principle, the non-resonant polarization term is also a 
sum of contribution from different wells, but since it does not depend on electric field, it is 
more convenient to keep it in the equation as a single term. 

Maxwell equation ([T]) together with polarization given by Eq. ([2]) and (jl]) describes a 
system of quantum well excitons interacting with common radiative field, E. In the ab- 
sence of polarization sources these equations have been intensively studied, and it is well 
known that they describe collective dynamics of radiatively coupled QW excitons^^'^Si'^E. The 




(3) 



UJmO -UJ -l-f 




(4) 
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main objective of the present section is to develop a general theoretical approach to solving 
Eq. ([1]), Eq. ([2]) and (jl]) in the presence of the sources of polarization of an arbitrary form. 
However, since in this paper the developed approach will be mainly applied to exciton lumi- 
nescence spectrum in the growth direction of the structure, we will restrict, for simplicity, 
our consideration only to s-polarized radiation. 

Using the translational invariance of the system in the x—y plane, we can present solutions 
of Eq. ([1]) in the form 

E(^,p) = e^'^''E(^,k), (5) 

where k is the in-plane wave vector. For a s-polarized wave the direction of k determines 
the direction of E(2;, k) as 

E(2;, k) = E{z, k)e,(k), e,(k) = e, x Gk (6) 

where e^, e^, and et are unit vectors describing directions of polarization, growth direction 
and the direction of the in-plane wave vector, respectively. In what follows we will omit the 
argument k when it is clear from the context that the value of the scalar amplitude is taken 
at a fixed value of the in-plane wave vector. 

Substituting Eq. ([5]) into Maxwell equation ([1]) and choosing s-polarized component of 
the field according to representation ([6]), we derive the following equation for the scalar 
amplitude of the field 

+n\z)E{z) = ^—F{z)- 



dz"^ c 



(7) 



dz' ^„,.{z')E{z') 



where n'^{z) = uj'^n'^{z) / — k"^, and S^, F{z) are the components of the two-dimensional 
Fourier transforms of the source terms Sm(p) and F(r) in the direction of e^: 



E^(k) =e,(k) • j d'^pHmip] 
F(k,z)=e,(k)- j d^p¥{p,z)e-''^f. 



(8) 



Equation ([7]) is the starting equation for the formalism developed below. The essential 
assumptions for this formalism are the non-local character of the exciton-light interaction and 
the possibility to separate in-plane coordinates. These assumptions are not too restrictive 
and, therefore, the formalism can be generalized for more complicated situations such as 



multi-resonance form of the susceptibility or the presence of the inhomogeneous broadening 
of excitons. The latter can be accounted for in the effective medium approximation by 
replacing the exciton susceptibility with its averaged over exciton frequencies version.— 



B. Transfer matrix approach to one-dimensional equations with sources 

The reduction of the initial problem to the one- dimensional equation (jTj) allows us to 
solve it using a powerful transfer-matrix technique. A convenient formulation of this ap 
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proach specifically adapted for the structures under consideration was developed in Ref. 
The presence of the source terms in Eq. ([7]), however, requires some modifications of that 
approach, and the adaptation of the transfer-matrix method to inhomogeneous integro- 
differential equations is one of the important technical results of this paper. 

Without any loss of generality we can consider a layer with the quantum well situated 
ai z = with the left and right boundaries at z_ and z+ respectively. Inside a single layer 
the summation over quantum wells in Eq. ([7]) as well as the well's index, can be dropped, 
and we can rewrite this equation in the form of a second order inhomogeneous differential 
equation, in which polarization terms appear as the right hand side inhomogeneity: 

^^ + k\z)E{z)=T{z). (9) 

A general solution of such an equation has the form^i 

E{z) = cih{z) + C2h2{z) + {G J^){z), (10) 

where /ii,2(-2) are a pair of linearly independent solutions of the homogeneous equation 

d^Ejz) 
dz^ 

and 



+ k'{z)E{z) = 0, (11) 



{G * J^) (z) = dz' G{z, z')T{z!) (12) 

Here G(2, z') describes the linear response of a passive (without exciton resonances) ID 
photonic crystal and can be expressed in terms of functions hx^2 as 

G(z, z') = Mz')h2iz) - hiz)hiz')] , (13) 

where Wronskian, = hidh2/dz — h2dhi/dz, does not depend on z. 
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We choose hi^2 as real valued solutions of the Cauchy problem for Eq. (fTTj) and use them 
to present the electric field at the left boundary of the elementary cell, z = Z-, a.s 



E{z. 



Cihi(Z-) + C2/i2(-2-)- 



(14) 



Combining Eq. ([7]) with Eqs. ffTOl) . f|T2|) . and (|T3|) we can derive the following expression 
for the value of the field at the right boundary of the elementary cell, z^: 



E{z+) =hi{z+) 
+h2{z+] 



Cl o fTTT- + X 3 I ClV^l + C2V?2 + 



2 I? A^. ,2, 



C2 + o /TT^ - X n I Clipl + C2V?2 + 



(15) 



where (/9i,2 and Fi^2 are the "projections" of the exciton state and the non-resonant field 
source onto the functions hi^2 

^= I dz<^[z)hi^2{z), 

(16) 



9^1,2 
Fl2 - 



1 



w, 



dz^{z)hi^2{z), 
dz F{z)hi^2{z)- 



In Eqs. (fT6l) the integrals are taken over the period of the structure (or over the elementary 
cell of the photonic crystal). The effective polarization source function S is the initial E 
modified by the field source function 

^ dz^{z){Gi.F){z). (17) 

In Eq. f[T5]) we also have introduced the modified exciton susceptibility 

X 



X 



Aujx/ " ' 



where 



Acu = ^^^/ dz^z){Gi.^)iz) 
c Jqw 



(18) 



(19) 



is the radiative correction to exciton susceptibility in the photonic crystal. 

Taking into account that the electric field at z = z^ can also be presented in the form 
of Eq. f|T^ with modified coefficients Ci^2, we can describe the evolution of the field upon 
propagation across the elementary cell of the structure as a change in these coefficients. 
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Using the solution (fTS!) the relation between the coefficients at different boundaries of the 
elementary cell can be found as 

(20) 

where the two dimensional vectors (ci,C2)(-2+) and (ci,C2)(-2_) represent the set of the re- 
spective coefficients, and Th is the transfer matrix describing their evolution across the 
elementary cell written in the basis of the linearly independent solutions hi 2'- 

r, = 1 + f V (21) 





i is the unit matrix. The contribution of the sources into the field is described by the second 
term in r.h.s. of Eq. ( l20l) 



(22) 



This is one of the main results of this Section. Its important feature is that the source 
contribution is independent of the state of the "incoming" field. In other words, the value of 
the field at the right boundary of the elementary cell is a superposition of a field propagated 
across the elementary cell from the left boundary (as if there were no sources at all) and 
the field generated by sources. This result is of a very general nature and can be applied to 
a variety of situations such as multipole exciton susceptibility, asymmetric quantum wells, 
incommensurability between periodicity of quantum well positions and modulation of the 
refractive index, etc. 

It might appear that Eq. ( l20l) violates the symmetry between the left and the right since 
the sources contribute only to the field at the right boundary of the elementary cell. This 
apparent asymmetry results from a fact that Eq. (!20|) presents a solution of the Cauchy 
problem, which is inherently asymmetric. The left-right symmetry should be expected only 
from a solution satisfying radiative boundary conditions and in the next section of the paper 
we demonstrate how to use Eq. fl2Ul) to find such a solution. 

The expressions presented in Eqs. ( 1211) and (l22l) can be greatly simplified in the case of 
symmetric quantum wells and the modulation of the dielectric function, which is consistent 
with this symmetry, i.e. n{zm + z) = n{zm — z), where Zm is the position of the center of 
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m-th quantum well. In this case the elementary cell of the structure can be chosen to have 
the explicit mirror symmetry with respect to its center (see Fig. [1]), and this is the case that 
we consider in what follows. 

Because of invariance of Eq. flTTl) with respect to mirror reflection, its solutions can be 
chosen to have a definite parity. Thus we can choose linearly independent solutions hi 2 to 
be either even or odd with respect to the center of the quantum well. For concreteness we 
choose /i2 to be the odd solution, which result in turning to zero and transfer matrix Tf 
taking the following much simpler form: 

In this case, the sources in Eq. f l22p can also be classified according to their symmetry. 
The amplitudes Fi 2 and (^1^2 have the meaning of projections onto the symmetric and 
antisymmetric solutions of Eq. ffTTj) . Respectively, only symmetric and antisymmetric parts 
of the sources contribute into these projections. The exciton radiative decay contributes 
only to Ac2 because the spatial distribution of the exciton polarization is determined by the 
exciton wave function, which is symmetric by the assumption. The non-resonant polarization 
F, in turn, does not have to have a definite symmetry and, therefore, generally contributes to 
both Aci and Ac2. In order to avoid any misunderstanding we have to emphasize, however, 
that the functions hi^2 do not represent the normal modes of an infinite photonic crystal. 
The latter are defined as solutions of an appropriate boundary problem and generally they 
do not have to be even or odd. 
As has been discussed in Ref. 24 



Eqs. ( 1201) and ( 1221) do not yet describe the propagation 
of the field across an entire elementary cell because they do not include the transfer across 
the interface between two adjacent elementary cells. The problem is that initial conditions 
for these solutions are defined at some point inside a given cell, and in order to use them to 
describe field in a different cell one has to introduce the shift of variables z — > z±nd, where d 
is the period of the structure, and n is the number of periods separating the two cells. After 
that one could express the functions with the shifted arguments as a linear combination of 
the original functions, but the most convenient way to describe the transition from one cell 
to another is to convert our transfer matrices to the basis of plane waves. In this basis the 
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FIG. 1: The periodic structure built of quantum wells (the shadowed rectangulars) and the barriers 
between them. Vertical dashed lines show the boundary of the elementary cell having the property 
of the mirror symmetry. The smooth line illustrates the modulation of the dielectric function in 
the structure. 

field and its derivative are represented as a superposition of waves propagating along ^-axis 

E = E+{z)e'''' + E_{z)e-''^' (24) 
dE/dz = iqE+{z)e'''' - iqE4z)e-''^' 

where q = k{z±), and can be naturally presented by a two-dimensional vector of the form: 

\E) = E+\+) + E_\-), (25) 



where 

are the basis vectors of the respective vector space. More detailed description of the plane 



1+) 



wave representation can be found in Ref. 
the amplitudes E± is written as 



24j . The relation between the coefficients ci 2 and 



^4- \ / Cl 

E^l \C2 



(26) 



where^^ 



M{z) 



hi(z)- 



iq 



ho(z) - 



iq 



(27) 



12 



Applying rule (126|) to Eq. (120!) we obtain 

\E){z+)=T\E){z.) + \vJ, (28) 

where T = M{z^)ThM~^{z^) is the transfer matrix through the entire period of the structure 
in the basis of plane waves. In the case of structures with the symmetrical elementary cell 
this transfer matrix can be presented as24 

af {af - a/)/2 
{af - fa)/2 af 

where 

a = 92, f = 91- iSg2, 
9*2, f = 9{ + iSg*2, 



T 



(29) 



(30) 



a 



and 



9i 



92 



1 



hi{z+) + 



KM 

iq 



'iqh2{z+) + h'2{z+)] . 



(31) 



The functions gi^2{'^)-, which are obviously not unique, are chosen to make clear the tran- 
sition to the limiting case of structures with spatially uniform refractive index. In this 
case choosing /?.i(0) = h^^Q) = 1 one has Wh = 1 and gi^2 = exp{iqz+). The function 
S{u!) = — 27rco'^xv^i/(9c^) introduced in Eq. (!30|) quantifies the interaction of the excitons 
with light. For the single-pole form of x(^) it has the form 

S = ^V— (32) 

uj — ujq — Auj + Z7 

where Tq = 2Tiauj'^(p\/ qc^ is the radiative decay rate. Because of the direct relation between 
the functions x ^i^nd S (they differ only by a factor slowly changing with frequency) we will 
for brevity refer to the function S{uj) as the exciton susceptibility. 
The source term (l22l) in the basis of plane waves takes the form 



\Vm) = M{Z+) 



Aci 

AC2 



-F2^- \a) + F,^- \s) ^ \s) , 



where 



la) = - 
' ' 2 





(33) 



(34) 



The index m in the notation \vm) reminds that all the relevant quantities can depend on 
the number of the well and should be taken for a particular well. 
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III. RADIATIVE BOUNDARY CONDITIONS AND THE FIELD EMITTED BY 
AN m-TH WELL 



Equation fl28l) expresses the field at the right boundary of the elementary cell in terms of 
the field given at the left boundary. Formally, it can be understood as the general solution 
of a Cauchy problem (the general solution of the homogeneous equation plus a particular 
solution of the inhomogeneous one). The amplitudes E± at the left boundary represent two 
independent parameters that can be chosen to satisfy any particular initial or boundary 
conditions. To represent a radiation coming out of the structure the field must satisfy 
radiative boundary conditions that require that outside the structure there must be only 
outgoing waves. Our objective now is, therefore, to find such Ej- that would satisfy this 
condition. To this end we consider an layer structure embedded into an environment with 
the refractive index Uout- Scattering of light by the interfaces between the terminal layers of 
the structure and the surrounding medium is described by the matrices 



Here zl^h are the coordinates of the left and the right ends of the structure, respectively. The 
angle of propagation is determined by ta.n6{z) = k/ n{z). The outg oing waves propagate at 
the angles following from Snell's law ?T,o„t sin 6'^^/? = n{zi^ji) sm6{zi^ji). 

We impose the radiative boundary conditions assuming first that the sources are localized 
only in the m-th layer. We require that in the half-spaces z < zl and z > zr, the field outside 
the structure would have the form of the wave propagating respectively to the left, and to 
the right. The former field can be described by a basis vector |— ) of the two-dimensional 
vector space introduced in Eq.( !25|) : E = E^^'^ |— ), and the later one is proportional to the 
other basis vector |+): E = E^^ |+). Using the results of the previous Section we can find 
the following relation between the fields outside the structure 




(35) 



where 



Pl,r 



Uout COS 9l,r - n{zL,R) COS 9{zl,r) 
Uout COS 9 L,R + n{zL,R) COS 9 {zl,r) ' 



(36) 



1+) = TRT{N,m + 1) \v^) + E'f^>TnT{N, 1)T-' |-) , 



(37) 



where 



T{N,m)=T^...T, 



m 



(38) 
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is the transfer matrix through the part of the structure obtained as a product of the transfer 
matrices through the individual layers. Eq. fl37|) is obtained by directly applying Eq. fl28|) 
to the field at the left end of the structure. This state is transferred through the entire 
structure in a usual way by a simple multiplication of the transfer matrices describing each 
period of the structure. This procedure results in the term proportional to the transfer 
matrix T{N, 1). Transfer across the luminescent layer results in an additional contribution 
as is given by the second term in Eq. fl28l) . After being emitted this field is then transferred 
across remaining N — m layers yielding the term proportional to T{N, m + 1). The matrices 
Tl^h take into account reflection of the radiation at the interface between the terminal layers 
of the structure, and the outside world. 

Multiplication of Eq. fl57|) from the left by (+| and ( — | gives the system of two inhomo- 
geneous equations with respect to E^^\ The solution of this system is 

rn(m) {-\TRT{N,m+ 1)\ 

^m/ 





HTpc 

\TLT-\m,l)\ 


Vm) 


(-1 


Tpc 


I-) ' 



(39) 



where Tpc = TjiT(N, 1)T£^ is the transfer matrix through the whole structure including 
the interfaces between the terminating layers and the surrounding medium. 

Equations ( l39l) can be used to derive expression for Green's function defined as a function 
relating the radiated field with the source. We can write the amplitudes of the waves outside 
in the form 



= Qi\m) (f, - S^a) ^ + Gt\m) (f, - t^^x) (40) 

where 

= TtN {-\TpT{N,m + l) \s,a) , g^^'''\m) = ±tjv {+\TLT-\m, 1) \s,a) . (41) 

Here we take into account the definition of the transmission coefficient through the whole 
structure t^r in terms of the transfer matrix ^at = (— | Tpc |— )~^- We would like to emphasize 
that these expressions are valid with the vectors |s) and \a) defined by Eqs. ( |34l) even for 
non-symmetrical structures where \s) and \a) do not imply symmetrical properties. 

Equation (l39l) can also be used to find the distribution of the field created by the source 
inside the structure. This can be achieved in two different ways. One can start, for instance. 
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with field Em^ |— ) on the left-hand side of the structure and propagate it across using 
transfer-matrices. When a luminescent layer is reached, Eq. |28] should be employed to 
describe transfer across it. Alternatively, one can propagate Em^ | — ) to find field in the 
elementary cells to the left of the luminescent layer, and propagate Em ^ |+) from the right 
to determine field in the cells to the right of it. 

We would also like to add that Eq. (HUj) can be used also for studying the directional and 
in-plane distributions of the field. For example, the standard problem of a point source can 
also be considered since the in-plane distribution of the source is taken into account by 2D 
Fourier transform ([8]). 

'as 
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Equation ( l40i) shows that the field outside is determined by two parameters founc 
convolutions of the sources with the functions hi^2{z)- As has been demonstrated in Ref. 
these functions provide complete description of photonic modes of a respective infinite sys- 
tem. From this perspective the result of Eq. fj40|) might seem expected. Indeed, it looks 
somewhat similar to a standard construction where the response is determined by a su- 
perposition of modes with amplitudes determined by the projections of the excitation onto 
the modes. This analogy, however, is misleading for the case under consideration. There 
are several important features distinguishing Eq. ( l40i) from the standard Green's function 
formalism. First, as has been noted, the functions hi 2{z) do not have to coincide with pho- 
tonic modes even locally (inside a particular elementary cell). Second, Eq. (HOj) is written 
for a finite structure when the applicability of the modes of an infinite structure can not be 
trivially justified. Third, Eq. ( I40l) does not require restrictive properties of the operators 
governing the propagation of light (say, hermicity) and, in particular, with slight modifica- 
tions remains valid in the presence of losses in the dielectric (Im[n(2;)] 7^ 0) while in this 
case even the notion of the projection has to be carefully examined. 

Below we will concentrate mostly on the case when the elements of the structure and 
the structure itself have the mirror symmetry. As a result, the index mismatch between the 
surrounding medium and the terminal layers is the same for both boundaries, so that one 
has Tl = Tr = Tp. It is interesting to note that in such structures the solutions given by 
Eq. f HOl) do not necessarily guarantee the symmetry of the radiation emitted to the left and 
to the right of the structure. For such symmetry to take place one has to prove that 
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This relation can indeed be proven but only in the case when both the structure and the 
sources are symmetrical with respect to the centers of the elementary cells. While the 
exciton related resonance source term does have the required symmetry, the non-resonance 
contribution is not necessarily symmetric if both F12 7^ (see Eq. (1551) and Eq. ([7j)). If this 
is indeed the case than radiation emitted to the right would not have the same characteristics 
as radiation emitted to the left. When F2 = 0, however, Eq. (H2!) can be proven with the 
help of relation \s) = —ax \s), where T is defined in Eq. f l29|) . and ax is the standard 
Pauli matrix. 

As follows from Eq. ( l40i) the field radiated due to the exciton recombination has, as 
expected, the resonant character. Thus, in resonant PCs the recombination determines 
the spectrum of the radiated field at frequencies close to t^o, while the no n- resonant sources 
specified by F{z) are responsible for the background component characterized by a relatively 
smooth frequency dependence. This contribution may become important farther away from 
the exciton frequency. For the problem of the exciton luminescence in the resonant PCs 
these non-resonant sources are not important and, therefore, below we assume that the only 
source of radiation is the exciton recombination and neglect the non-resonant contribution. 
The expression for the radiated field essentially simplifies in this case and takes the form 

E^r^ = -^±(m)S^^, (43) 

where we explicitly show the dependence of the source on the number of the layer and 

(m) = ±tM (± \T^im)T,\s) . (44) 

Here we have taken into account that = when F = and have dropped the super- 
script (s) since the only relevant Green's function in this case is g±\m). We have also intro- 
duced partial transfer matrices T^{m) = TpT^^{m, 1)T^^ and T+(m) = TpT{N,m + 1)T^^, 
which have the property TZ^{m)T+{m) = Tpc- 
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IV. THE LUMINESCENCE SPECTRUM OF RESONANT PHOTONIC CRYS- 
TALS 



A. Quasi-classical approach to luminescence 

In this section we will apply general results of previous sections to the problem of the lu- 
minescence spectrum of multiple-quantum-well based resonant photonic crystals. Lumines- 
cence is one of the manifestations of spontaneous emission, and as such is a purely quantum 
electro dynamic phenomenon. At the same time, the non-coherent radiation produced due 
to luminescence in many situations can still be described as classical electromagnetic field 
with randomly changing amplitude and phase. The average value of such a field is equal 
to zero, while its variance, defined as (£'^(r)), is identified with intensity, where angular 
brackets (■ ■ ■ ) symbolize averaging over appropriately defined distribution function. Such a 
quasi-classical approach to spontaneous emission has been used previously in a great number 



of different situations. Just as a small sample of papers dealing with applications o: 
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classical 
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electrodynamics to the problem of spontaneous emission we can refer to Ref. 

One of the methods of reproduction of non-coherent classical field is a Langevin-like ap- 
proach, in which the polarization sources appearing in macroscopical Maxwell equations 
([1]) are considered as random functions of time and coordinates whose statistical charac- 
teristics should be determined either from experiment or from fully quantum microscopical 
theory. In what follows we will neglect the non-resonant contribution to polarization F, and 
characterize statistical properties of the resonant source function Sm(p, by a correlation 
function 

(Si,^(p, t)S,, i(p', t')) = SijSmiKY: (p - p', t - t') , (45) 

where (■ ■ ■ ) signifies statistical averaging over various realizations of the non-coherent exciton 
polarization, indexes i and j designate Cartesian coordinates in the {x, ?/)-plane, and m 
and / are well numbers. Equation ( H6l) implies that the fluctuations of the non-coherent 
exciton polarization are (i) statistically uniform in time and space, (ii) direction of the non- 
coherent polarization is distributed isotropically in the plane of the structure, with various 
components of the polarization vector independent of each other, and (iii) source functions 
in different wells do not correlate with each other. 

This correlation function can only be found from a microscopical theory of electron-hole 
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relaxation processes. An example of such a theory, which provide a more quantitative jus- 
tification for our phenomenological approach and shows relation of this correlation function 
to microscopic characteristics of quantum wells, is presented in Appendix. These calcula- 
tions demonstrate that all three assumptions regarding properties of the source correlations 
are justified. The most important of them is the assumption of independence of the source 
functions in different wells. We would like to emphasize here that this assumption does not 
mean that we neglect radiative coupling between wells. As it is explained in Appendix the 
source term Sm(p, is determined by excitons populating non-radiative states so that the 
electromagnetic field associated with them decays exponentially outside of the wells and 
does not contribute to the radiative coupling. This coupling is described by the common 
radiative field E which is created by and act on all quantum wells in the system as it was 
explained in previous section [III It is interesting to note that Maxwell equation ([1]) describes 
this coupling even though the average value of the non-coherent field is equal to zero. The 
correlation function for the fourier transformed source function T,m{uJ, k) is given by spectral 
density k) defined as 

(S„(k, uj)Ei{k', Lo')) = E^oo, k)5{uj - ^')^(k - k')<5„z (46) 

which is a Fourier transform of the time-position correlation function given in Eq. fHSj) . 
The field created by such source is characterized by a spectral intensity X(k, uj) defined 

as 

(£;(ki, cui)E(k2, u^)) = T(ki, u;i)5(ki - k2)5(u;i - u^). (47) 

Applying Eq. (H6!) to this equation we find the spectral intensity of radiation emitted by the 
entire structure in the form 

T±iuj,k) = 4 ^Em{uj,k)\g±{m;u,k)\'^ 

m 

which implies that the field emitted by different wells adds in a non-coherent way. This 
general expression allows analyzing both the frequency and the directional dependence of 
the luminescence spectrum. In this work we restrict our consideration to the waves emitted 
along the growth direction of the structure (i.e. k = 0). The directional distribution of the 
radiation will be studied elsewhere. 

Equation (H8|l shows that the form of the luminescence spectrum is determined by several 
factors with different frequency dependencies. The exciton susceptibility, S{uj) for instance, 
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qSm{uJ, k) 



(48) 



strongly reduces the luminescence far away from the exciton frequency ooq- The spectral 
density H(ti;) is expected to show a weak frequency dependence at the scale of the width 
of the polariton stop-band. The factor \Q {m; u , k)]"^ , according to Eq. (H^ . is the product 
of two terms. One is the transmission coefficient tjv? which may have strong frequency 
dependence following the singularities at the eigenfrequencies of the quasi-modes of the 
structure. These singularities determine the fine structure of the luminescence spectrum. 
The second term is responsible for the variations in the luminescence intensity at a much 
larger scale. Equation (148|) can be used to analyze spectra of luminescence of various types 
of MQW structures. Several examples are presented in the subsequent parts of this section. 



B. The luminescence spectrum of finite periodic structures 

1. General expression for intensity of emission 

We first apply our general results to the case structures built of identical layers, so that 
V^mi = V'l and Sm = S, i.e. all quantum wells are characterized by the same exciton frequency 
uq. We would like to note that in order to have the structure with the mirror symmetry 
it must contain an integer number of the elementary cells. This means, in particular, that 
the terminating layers are half-barriers. As has been noticed, in this case the luminescence 
spectrum is the same at the both sides of the structure. For the concreteness, we will 
consider the field radiated to the right (i.e. E^). 

In the symmetric case we can perform the summation over the quantum wells using the 
fact that all partial transfer matrices posses the mirror symmetry and can be presented in 
the forn>2i 

(cos 9 — i sin 9 cosh B —i sin 9 sinh 3 \ 
h (49) 
i sin 9 sinh /3 cos 9 + i sin 9 cosh /3 / 

where 

cos9 = ]-{af + af), tanh /? = ^^^^^^ — (50) 
2 af -af 

The parameter 9 determines the polariton spectrum of an infinite structure and is defined 
as 6' = Kd, where K is the polariton Bloch wave-number. 

The representation fj49l) is convenient because all T^{m) are characterized by the same 
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/5, while the spectral parameter of T^{m) is merely —mO. Thus, T^{m) can be written as 



Ui(3/2)\+) {+\U-\P/2) 



+ e 



-imd 



U{P/2)\-) {-\U~\P/2), 



(51) 



where U{l3/2) = TpTuifi /2), and Th is a matrix describing a hyperbolic rotation with a 
dilation 

„ , cosh [3 — sinh /? 



sinh 13 cosh /3 



(52) 



Matrix Tp takes into account reflection and transmission at the external boundaries of the 



system. 

Using Eq. (!5T!) one can flnd 

v^,^ . .,9. sinhNO" 



\tN\- 



sinh ^" 



|^|2g-(iV+l)e" ^ 




sinA^^' r 






sin 9' . 



(53) 

where we have introduced real and imaginary parts of the dimensionless Bloch number, 9: 
9 = 9' + i9" and parameters 



A 



B 



1+P 
1 



{92 - pgl) cosh^(/3/2) - -{gl - p^s) sinh/3 
(^2 - pgl) sinh2(/3/2) - ^(^; - p^2) sinh/3 



(54) 



For the purposes of numerical calculations instead of direct calculations of the parameter 
(3 it is more convenient to multiply both parts of Eqs. by sin 9 and, then, use Eq. (149|) 
to establish direct relation of corresponding terms with the elements of the transfer matrix 
through the period of the structure. 

While Eq. (1531) describes luminescence of a periodic MQW structure with an arbitrary 
period, we shall focus our attention to the most interesting cases of Bragg and near-Bragg 
structures,— in which effects of periodic modulation of the refractive index and light- 
exciton coupling are most pronounced. As we already mentioned in Introduction the period, 
c?, of such structures satisfles a special resonance condition ojq = u!B{d), where the exact 
value of the resonant frequency ub in systems with periodically modulated refractive index 
depends not only on the period of the structure, but also on details of the modulation.—"^ 
For concreteness we will assume that the dielectric function reaches its maximum value at 
the quantum well and monotonously decreases towards the boundaries of the elementary 
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FIG. 2: (color online) Evolution of polariton band boundaries with detuning of the structure with 
piecewise modulation of the refractive index from exact Bragg condition. The horizontal axis 
describes the detuning in terms of parameter to — ujq, where w = TTc/{nwdw + ribdh), n^, nh, 
dh and refractive indexes and thicknesses of well and barrier layers respectively. The vertical axis 
shows a relative difference between frequency uj, and exciton frequency ujq. 

cell. In this case the Bragg resonance takes place when the exciton frequency coincides 
with the high-frequency boundary of the photonic band gap, Since the details of the 
emission spectrum are determined to a large extent by the electromagnetic band structure 
of the systems under consideration, it is useful to remind the main features of this structure, 
which was analyzed in details in a number of papers.— i^^iSiiS^iS^ Fig. [2] shows the dependence 
of band boundaries of MQW structure upon its period, where shaded regions correspond to 
polariton stop-bands. One can see that at a certain value of Co = 7rc/{nwdyj + ribdh), where 
riw, fib, dyj, db and refractive indexes and thicknesses of well and barrier layers respectively, 
two stop-bands connect at the exciton frequency ujq forming a single wide band-gap. This 
is the point of the Bragg resonance, when exciton frequency falls inside a stop-band of 
the spectrum, whose width can be much larger than the width of the exciton resonance. 
(The second occurrence of a single band situation at larger values of lj results from collapse 
of one of the gaps, and is a result of random degeneracy between two exciton polariton 
branches.) Spectrum of structures only slightly detuned from the Bragg condition (we will 
use term quasi-Bragg for such structures), is characterized by emergence of a propagating 
band between the two stop-band. The exciton frequency in this case belongs to the boundary 
of the propagating band, which is situated asymmetrically with respect to the outer band 
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boundaries: for negative detunings ujq is closer to the upper boundary, while for positive 
detuning the lower polariton branch eventually moves closer to it. 

Equation( [55]) demonstrates how the polariton band structure affects the luminescence 
of the system under consideration. One can see from this equation that the structure of the 
spectrum is characterized by two scales of frequencies. On a smaller scale the modulations of 
the intensity of emission are determined by |tivp term in Eq. (1531) . whose maxima correspond 
to the real parts of the polariton eigenfrequencies. The modulation of intensity on this scale 
depends on the number of periods in the structure and occurs over frequency intervals of 
the order of Vg/{dN), where Vg is the group velocity of the polariton excitations and d is 
the period of the structure. The polariton band structure affects the luminescence on a 
much larger spectral scale through the combination of the exciton susceptibility S{uj) and 
the imaginary part of the polariton Bloch number presented by parameter 6" . 

Homogeneous and inhomogeneous broadenings significantly effect the short-scale modu- 
lations of the luminescence allowing for their observation only in high quality samples at 
very low temperatures. At higher temperatures the long scale variations of intensity, which 
depend significantly on relations between ujq and ub, become predominant. In the case of 
Bragg structures, when uq is very close to ub, Eq. fl53|l predicts that the luminescence spec- 
trum is mostly concentrated outside of the polariton stop-band near the edges of the bands 
of the exciton polaritons. Indeed, at frequencies inside the forbidden gap the contribution 
to I{uj) of the exponentially large terms in the right-hand-side of Eq. (!53l) is canceled by 
the exponentially small transmission at these frequencies. As a result, only wells within 
the attenuation length from the boundaries contribute to the radiated field. Besides at fre- 
quencies far away from the ujq the luminescence is subdued by the smallness of the exciton 
susceptibility. These qualitative arguments can be supported by direct calculation of the 
emission intensity in the neighborhood of the band edges, where Eq. (1531) can be simplified. 
In this spectral region we can represent the spectral parameter a.s 6 = tt + ie and assume 
that e is sufficiently small, so that A^|e| ^ 1. Obviously, this approximation covers a sub- 
stantial interval of frequencies only for not very long structures, but it is quite sufficient for 
a qualitative analysis of realistic structures. 

Expending Eq. (l48l) in a power series with respect to the small parameter eN, we can 
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FIG. 3: (color online) The luminescence spectrum in comparison with a polariton band structure 
is shown for quasi-Bragg structures consisting of 60 layers near the boundary of the first Brillouin 
zone. Smooth filling corresponds to the luminescence calculated using H = 1. The lines are the 
contour plot of 9" , where the outermost curve correspond to its smallest value. Frequency changes 
along the vertical axis, while, the horizontal one presents detuning from the Bragg resonance, (a) 
Pure MQW structure with parameters typical for Ala;Gai_i.As/GaAs structures: Fq = 75 //eV, 
loq = 1.489 eV, 7 = 300 /xeV. (b) An example of MQW based photonic crystal. The exciton 
related parameters are the same as in (a). The modulation of the index of refraction is taken to 
be n{z) = 3.4 + 0.1 cos20(7rz/2d). 



approximate it by the following simple expression: 

S{oo) ' hlq' 



(55) 



An important result immediately demonstrated by this equation is a linear increase of the 
intensity of the emission with the number of quantum wells. This is an expected behavior 
because of the transparency of the structure at these frequencies and the independence of the 
contributions of different wells to the emitted light. Another important conclusion following 
from Eq. fl55l) is a relative weakness of the luminescence of the Bragg structures. The cause 
of the decrease in the emission is related to the presence of a broad polariton stop-band in 



such structures whose width, A, given by expression A = ^/2r\pj^jTi , is much larger than the 
width of the exciton susceptibihty Siuj) determined by non-radiative decay, 7. In the interior 
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of the stop-band the luminescence is suppressed by small transmissivity of the structure in 
the vicinity of the exciton resonance, while at the edges of the stop-band it is reduced by 
the factor of Fq /cuq ^ 1 because of the separation of the band boundaries from the exciton 
frequency. The found decrease in luminescence for Bragg structures is equivalent to the so 
called sub-radiance effect obtained in Ref. 
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on the basis on fully quantum calculations. 



We demonstrate here that this effect has a purely classical origin and is caused by formation 
of polariton stop-band in Bragg MQW systems. 

Detuning from the Bragg resonance opens up a transparency window inside the stop-band 
with exciton frequency coinciding with one of the band boundaries (Figl2]). On the base of 
the same arguments as above one can expect that the emission spectrum is characterized 
by two maxima: the stronger one in the vicinity of the boundary of the propagating band 
adjacent to ujo, and the second, weaker, maximum at the boundary of the outer polariton 
band, which is closer to the exciton frequency. The second outer boundary of the polariton 
band is so remote from the exciton frequency that its contribution to emission can be 
neglected. Numerical calculations carried out with the exact form of I{uj) confirm these 
conclusions. The results of these calculations are presented in Fig. [3l where we show the 
dependence of the intensity upon frequency and the period of the structure. The intensity is 
shown by the shading on the graphs — the darker shading corresponds to higher emission. 
In order to facilitate a better understanding of the role of the refractive index contrast 
we simulated two types of structures: one with a realistic changes in the refractive index 
between wells and barriers, and the other, in which refractive index was assumed constant 
throughout a structure. The latter structures are often called optic lattices because all the 
modifications in their optical properties come from the radiative coupling between quantum 
well excitons. It is interesting to see a significant difference between the luminescence spectra 
of MQW optical lattices and MQW based photonic crystals. The latter is asymmetric with 
respect to the point of the Bragg resonance, while the former shows complete symmetry. 
This feature is clearly related to the asymmetrical structure of the polariton band gap in 
structures with modulated refraction index.-^ In order to emphasize the relationship between 
the luminescence spectrum and polariton band structure, the spectrum in this figure is 
presented together with the polariton stop band. The latter is shown with the help of 
level curves of the imaginary part of the dimensionless Bloch vector 6" . In an ideal system 
without any broadenings, 6" would have been zero everywhere outside of the stop-band. In 
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(a) (b) 

FIG. 4: The fine structure of the luminescence spectrum. The parameters of the structures are the 
same as in Fig. [3] except 7 = 30 /ueV. For better visibility the shadow intensity is chosen according 
to the log-scale, (a) The MQW structure with a homogeneous dielectric function, (b) The MQW 
based photonic crystal. 

real systems, of course, 6" is not zero everywhere because of the exciton broadening. This 
makes the notion of the stop-band not very well defined, and, in particular, the edges of 
the gap can not be determined unambiguously. However, at the frequency, which would 
correspond to the band edge in a system without broadening, the imaginary part of the 
polariton Bloch wave-number drastically increases. This increase can be traced on the level 
curves of 6" in Fig. [3], where outer curves correspond to the smallest value of 6". It is seen 
that the maxima of the luminescence spectrum approximately follow these lines when the 
relation between the the exciton frequency and the period of the structure changes. The 
exact position of the maxima is determined by an interplay between a smaller value of 6" 
(and, hence, a higher transmission) and a smaller distance from the exciton frequency (a 
higher value of S{uj)). Comparing the spectrum shown in Fig. [3] with the band structure 
shown in Fig. [2]one has to notice the different frequency scales of these figures. The frequency 
region covered in Fig. [3] includes only the small transparency window around the exciton 
frequency and only the closest to it outer polariton band. 

For sufficiently smaller value of the exciton broadening the fine structure becomes clearly 
visible as is seen in Fig. HI As we mentioned above the maxima of the luminescence forming 
this fine structure result from the periodic dependence of the transmission on frequency. 
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FIG. 5: (color online) Fine structure of the PL spectra for the 100 wells MQW structures with 
different periods. Curves 1, 2 and 3 are calculated for the period respectively smaller, equal, and 
larger than the Bragg resonance value. The indexes of refraction of wells and barriers in this 
calculation were assumed equal to each other. 

These maxima appear as the characteristic scars on the spectrum presented in this figure. 
More clear representation of these features of the luminescence spectrum can be given by 
direct plotting of the intensity as function of frequency for different values of the de-tuning 
of the structure from the Bragg resonance, as shown in Fig. 0, which was obtained neglecting 
the modulation of the refractive index. 

2. Comparison with experiment and the role of inhomogeneous broadening 

Comparing our calculations with experimentally observed spectra,^^'^^ one should take a 
few considerations into account. First of all, the direct quantitative comparison is rather 
difficult because the experimental spectra are infiuenced by details of the entire experimental 
sample, and not just by its MQW part. For instance, the details of the cladding layer can 
significantly influence the observed luminescence spectrum. In order to illustrate this point 
we used the general formulas derived in the paper to calculate the emission intensity in the 
presence of the cladding layer. The results of these calculations are shown in Fig. El where 
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FIG. 6: (color online)The same spectra as in Fig. [5] but with cladding layer of thickness dc = 
db + dw/2 

one can notice significant changes introduced by the cladding layer to the spectra. Second, in 
photoluminescence experiments with long MQW structures the intensity of pump radiation 
is not uniform along the structure, which results in different source functions for different 
wells. This circumstance also affects observed spectra as can be demonstrated by direct 
computations using general formulas obtained in this Section. To this end we assumed that 
the source function, Em{uj), which appears in Eq. PHI) can be presented as an exponentially 
decreasing function of the well number, m: Sm(ti^) oc exp{—amd), where parameter a 
represents an inverse attenuation length of the pump. Using this representation for the 
source function in Eq. (HHl) we numerically calculated emission intensity with different values 
of parameter a. The results of these calculations are shown in Fig. [TJ where luminescence 
spectra with a = and a = 0.2 are compared. This figure clearly demonstrates that 
inhomogeneity of the source function can have a significant impact on the observed spectra. 
Having in mind mentioned circumstances, we will not attempt to quantitatively reproduce 

eatures^jwhich most likely 
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and 
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one can notice 



experimental spectra, focusing instead on the most significant 
have intrinsic origin. Comparing experimental results of Refs. 
that despite of quantitative difference between these two experimental spectra, they share 
one common feature, which is, at the same time is in a striking contrast with results of 
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FIG. 7: (color online) The luminescence spectrum of the almost Bragg 100 layers structure with 
exponentially decaying source term characterized by the decay rate a: curve 1 corresponds to 
a = 0, curve 2 to a = 0.1 and curve 3 represents a = 0.2 

our calculations. According to our predictions, the luminescence must be most intense in 
the vicinity of the exciton frequency, while the experiments show that out of two most 
pronounced maxima of the emission, the one, which is farther away from ujq is brighter. 

One probable reason for this discrepancy is the inhomogeneous broadening of excitons, 
which has not yet been taken into account in our calculations. In this work we include effects 
due to the inhomogeneous broadening into consideration using a simple model of effective 
medium.—"^ Within this model one neglects spatial dispersion of excitons and assumes 
that inhomogeneous broadening is caused by spatial fluctuations of exciton frequency ujq. 
It is further assumed that these fluctuations can approximately be taken into account by 
replacing exciton susceptibility, Eq. ([3]), in all relevant equations, with its average value 



where p{z) is a distribution function of exciton frequencies. In the case of a not very strong 




(56) 



29 



((Oq-COb)/© 



-0.02 -0.01 0.00 0.01 0.02 



FIG. 8: The same spectra as in Fig. [3^, but with inhomogeneous broadening taken into account 
within the effective medium approximation. The parameter of inhomogeneous broadening was 
chosen to be equal to a = 200 jieV. 



inhomogeneous broadening this function can be approximated by a Gaussian: 

\2" 



1 



exp 



[Z - UJq) 



(57) 



where ujq is an average exciton frequency, and its r.m.s. value, a determines the width of the 



distribution. This model of the inhomogeneous broadening was first introduced in Ref. 
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on heuristic basis for calculations of refiection and transmission spectra of MQW structures. 



and later more rigorously justified in Ref. 
Ref, 
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While derivation of this model carried out in 



cannot be directly applied to the luminescence problem, we still believe on the ground 



40, that the effective medium 



of physical arguments similar to those put forward in Ref. 
approach can give qualitatively accurate description of the role of inhomogeneous broadening 
in luminescence spectra, at least, for emission in directions close to normal. Indeed, light 
emitted by quantum well excitons in the close to normal direction leaves the quantum 
well without moving significantly in the in-plane direction, and thus without experiencing 
significant scattering due to in-plane disorder. Accordingly, the main effects of the in-plane 
disorder in this case is that excitons localized in different regions of the sample emit light 
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FIG. 9: Cross-sections of the intensity profile shown in Fig. [8] for three different values of detuning 
from exact Bragg resonance. 

at different frequencies. A probe with a sufficiently large aperture (a typical situation in 
luminescence experiments unless one deals with microluminescent spectra) would collect 
light emitted by all these excitons, effectively averaging out their susceptibility. 

The result of numerical computation of emission spectra with Eq. (!56|) for exciton sus- 
ceptibility are presented in Figs. ([8]) and (Q. The first of these figures show that taking into 
account the inhomogeneous broadening resulted in spectral redistribution of the emission 
intensity from peaks closer to the central exciton frequency to those that are farther away 
from it. The latter are now brighter than the former in a qualitative agreement with exper- 
imental spectra. This point is demonstrated even more clear in the second of these figures, 
which shows emission spectra for three values of detuning from Bragg resonance. Qualita- 
tively this affect can be understood by noticing that by averaging the exciton susceptibility 
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we essentially smoothed its resonance dependence on the frequency reducing, therefore, ef- 
fect of decreasing susceptibility on the emission intensity. In this situation the intensities of 
luminescence peaks are determined by an interplay between affects due susceptibility and 
transmissivity of the structure. These calculations show that inhomogeneous broadening 
can be in principle responsible for observed luminescent spectra, while it is clear that a 
quantitative agreement with experiment would require more rigorous treatment of inhomo- 
geneous broadening as well as taking into account such effects as inhomogeneity of pump 
and cladding layers. 

C. The luminescence spectrum of structures with defects 

One of the main reasons for the interest to resonance photonic crystals derives from 
possibilities to manipulate their optical properties through modification of their structure. 
One of the possible approaches includes intentional violation of the periodicity of the PCs 
by introducing one or several defects. In one-dimensional structures such defects are lay- 
ers with different characteristics. Depending upon which parameters of the defect layer 
are modified one can have a variety of defect structures. Effects of such defects on re- 
flection/transmission properties have been studied for both regular— >^!^!^ and resonant 
photonic crystals.— 1^'^'^'^'^'^ In the particular case of Bragg MQW structures it was 
demonstrated, for instance, that by using different types of defects one can engineer struc- 
tures with a wide variety of optical spectra.^^ It is clear that defects will also substantially 
affect luminescence properties of the Bragg structures. Despite the obvious interest of this 
issue for applications, it has not yet been addressed, and in this paper we present the initial 
analysis of this problem for one particular type of the defect structure. 

The role of the defects on the luminescence spectrum of Bragg MQW structures is two- 
fold. Firstly, the defects affect the emission of the regular part of the structure caused by the 
modification of the transmission spectrum, tjy, of the structure. Secondly, the defect layers, 
depending on their structure, can contribute their own luminescence to the total spectrum. 
It should be noted, however, that the second contribution is expected to be small because of 
the much smaller number of the defect layers compared to the total number of the periods 
in the structure. 

In this paper we will illustrate a possibility to modify luminescence of Bragg MQW 
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structures with the help of structure manipulation by considering one particular case of a 
defect structure. We will consider a (2A^ + l)-layer multiple quantum well structure, in 
which the first and the last layers have a fixed width w, while (A^ + l)-th layer at the 
middle does not contain a quantum well and has a width d. To simplify our analysis we will 
neglect the refractive index contrast between well and barriers in the structure. This type 
of defect can be described as a cavity, in which parts of the structure to the right and to 
the left of the defect layer are identified as mirrors. Transmission and reflection properties 
of this structure in the region of the stop-band have been studied in Ref. |45| in the limit of 
very long structures. Here we will consider more realistic case of relatively short structures, 
and will analyze the manifestations of this defect in luminescence. 

The transmission properties of the structure are described by the transfer matrix 

T = TjTrfTj, (58) 

where and are transfer-matrices through layers with and without quantum wells, 
respectively, 

(j<l>^(l_iS) -iS \ (e^'^'^ \ 

iS e-"t'-'{l + iS) j I e-"l''i I 

with (f)y^ = kw and (pd = kd. 

The luminescence spectrum can be calculated using general expression (HHl) . As has 
been noted, the strongest effect on the luminescence spectrum can be expected due to the 
modification of the transmission by the defect layer. We consider this effect for the case 
when it is most pronounced, i.e. when MQW parts of the structure satisfy the Bragg 
condition, 0^„(co'o) = tt. In A^-layer MQW structures without the defect the transmission 
in this case has a deep with the width ~ NTq centered at the exciton frequency, c^o- The 
resonant tunneling induced by the defect mode results in the appearance of the resonant 
transparency that, in turn, leads to resonant grow of the luminescence at the respective 
frequencies. In order to provide a qualitative description of this effect it is more convenient 
to work with reflection r = — T21/T22 and consider its resonance drop caused by the defect 
layer. Using the representation (H9|) for the transfer matrices surrounding the defect layer, 
one can find for not too long structures 

T21 cx cos (pd + N sin 0rf(sin (p^ — S cos (p^). (60) 
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FIG. 10: (color online) (a). The luminescence spectrum of the Bragg 40 layers structure with the 
defect layer inserted at the middle of the structure. The intensity of the emitted field is plotted 
as a function of frequency (the vertical axis) for structures with the defect layer characterized 
by different LOd (the horizontal axis). The solid line depicts the position of the maximum of the 
transmission as obtained in Eq. (j6ip . The parameters of the structure are the same as in Fig. [3^. 
(b). Explicit form of the luminescence spectrum of the defect structure for several values of the 
defect resonant frequency uid- 

Taking into account that T22 does not have a significant frequency dependence, we can 
assume that the minimum of the reflection, uOmim occurs at the same frequency as the 
minimum of T21. Neglecting the homogeneous width of the exciton resonance compared to 
the width of the stop-band Ap, we can present an equation determining Umin in the form: 

1 arctan I , (bl j 



where ojd is defined by (pdi^d) = tt. 

As follows from this consideration, the defect layer inserted into a MQW structure leads 
to the appearance of a transparency resonance with the width ~ 7. The position of the 
maximum of the emission is determined by the interplay between the maximum of the 
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transmission and the decay of the exciton luminescence at frequencies away from cuq. The 
maximum of the field intensity, therefore, is expected at frequencies shifted from Umin to- 
wards Uq. 

In Figure [TD] we plot the luminescence spectrum for the structures with different width 
of the defect layer. The spectrum is obtained by performing direct summation in Eq. 
using the invariant embedding method^ for finding the partial transfer matrices. As one can 
see, the numerical calculations confirm the simple analysis provided above and demonstrate 
a sharp rise in the luminescence due to the defect state. 



D. Absorption and luminescence spectra 

According to the Kirchhoff 's law all emitting system in thermodynamic equilibrium must 
demonstrate an universal relationship between their absorption and emission spectra. Such a 
relationship exists even in non-equilibrium but stationary situations such as luminescence in 
the regime of steady state excitation. To establish such a relationship for a particular system, 
however, is not always a straightforward task. For instance, authors of Ref. [sj managed to 
establish such a relationship for quantum well excitons only as an approximation. At the 
same time, this relationship might be rather important because it could allow to relate to 
each other various microscopic characteristics of a system under consideration independently 
of particular microscopic model used for their calculation. Here we will use the approach 
to description of resonant photonic crystals developed in the present paper as well as in 



Refs. 
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l37l in order to derive an exact relation between absorption and emission spectra of 
resonant photonic crystals. 

A traditional definition of absorbance, A{uj), of an open one-dimensional dielectric struc- 
ture has the form 

A{uj) = 1 - T{uj) - R{uj) (62) 

where the last two terms represent transmission and reflection coefficients respectively. This 
form can be effectively used to carry out a thermodynamical derivation of the relation 
between emission and absorption specifically tailored for one- dimensional systems. Let's 
assume that our one-dimensional structure is bounded by vacuum on its left-hand side and 
by a homogeneous dielectric medium with refractive index rit on its right-hand side. The 
derivation is based on the statement that in equilibrium total photon fiux out of the system 
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must be equal to zero. This total flux includes the flux of emitted photons, $em, and the 
flux of the incident and transmitted/reflected photons, $i, $j, and $r, respectively. $em 
can be calculated as 

= h{u)du^^, (63) 

where /q is the equilibrium luminescence intensity, and q^^y represent conserving in-plane 
components of the photons' wave vector. The incoming and reflected fluxes in vacuum can 
be presented as 

dqydq^ 



= hooN'^^ioo) [1 - R{uj)] doo^^ , 

where Npf^{uj) is the equilibrium photon occupation number. The last contribution to the 
flux of photons in vacuum comes from the photons transmitted from the medium on the 
right-hand side of the structure. This flux can be written down as 

= ^Kd^)nco)^^0^ = hu^Nl,{uj)f{u^)du^ (65) 

where T stands for the transmission coefficient of light incident on the system from the 
right-hand side, and we took into account the refractive index nf, of the medium on the 
right. Combining Eq. fl^^ with Eq. fl64p . fl5S]) and with the requirement of the zero total flux 
we can write down 

hujNly^{uj) [l - R{uj) - f[iA = Jo(cu) (66) 

Taking into account Eq. fl62l) and the fact that due to the time reversal symmetry transmission 
coefficients, T for the wave incident from right is equal to the one describing waves incident 
from left we obtain a final form of the Kirchhoff's law for one-dimensional layered structures 

loiiu) = huN'^,iuj)Aiuj) (67) 

Equilibrium photon distribution function at low temperatures ksT -C hcu can be approxi- 
mated by exp {—hu/kBT). 

The derivation of Eq. fl57|) is based on two main assumptions: zero photon flux, and 
time reversal symmetry, both of which are valid for any kind of steady state, not necessarily 
equilibrium, situations. In non-equilibrium cases, however, the concept of photon distri- 
bution function, Nph{uj), appearing in Eq. (l67|l . becomes ambiguous, and generalization of 
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this equation to these situations is not straightforward. Nevertheless, under certain circum- 
stances, a relation similar to Eg. (1671) can be obtained for such non-equilibrium phenomena 
as luminescence under steady state excitation. Results of numerous studies of exciton photo- 



luminescence in QW (see, for instance, Ref j27l]5ll ) indicate that in the case of non-resonant 



photo-excitation of luminescence at the moderately low temperatures where the main con- 
tribution to luminescence already comes from free excitons, the following kinetic of lumi- 
nescence can be assumed. Originally excited electron-hole pairs first relax through phonon- 
assisted processes to high energy states with in-plane wave numbers k corresponding to 
non-radiative " dark" excitons^i. These states live long enough to come in quasi-equilibrium 
with the crystal lattice, so that they can be characterized by a Boltzmann distribution 
fsiE) oc exp [—{E — fi)/kBT]. The dark exciton act as a source of exciton luminescence, 
and, according to calculations presented in Appendix, the photoluminescence intensity I{uj) 
is a linear functional of the exciton distribution function, fsiE). Term explfi/ksT) in 
the Boltzmann distribution can be factored out, and the remaining expression reproduces 
equilibrium emission intensity. Thus for the intensity of the luminescence we can write 
I{u) = exp {fi/kBT)Io{u), from which it follows that /(a;) and A{uj) are related to each 
other as 

I{uj) = hiuexp {^i/kBT)Nl^{uj)A{uo) (68) 

Experimental verification of the relation given by Eq. (|68l) can provide useful qualitative 
information about the distribution of excitons and their kinetics either this relation is con- 
firmed or not. Eq. fl68|) can be applied to emission and absorption of each QW constituting 
the structure, and if the distribution functions of excitons are identical in each well, to the 
entire MQW structure as well. 

The established relation between emission and absorption is based on thermodynamical 
arguments and depends little on the details of the structures under consideration. It appears 
useful to compare this relation with the one derived on the basis of solution of Maxwell 
equations for the particular model of MQW structure considered in this paper. To this end 
it is more convenient to use an alternative expression for absorption, which follows directly 
from the definition of Poynting vector and energy conservation: 

where S is the Poynting vectors of outgoing radiation respectively, and L^, Ly are the 
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magnitude of the Poynting vector of the incoming radiation and the transverse dimensions of 
the sample respectively; the integral is taken over a surface enclosing the entire sample. The 
particular convenience of Eq. (1691) for ID structures stems from the fact that this expression 
allows for expressing the absorption in terms of the fields at the boundaries of the structure. 
Assuming that the incoming radiation impinges on the structure from the left, where the 
space is filled with the medium with refractive index n^, one can obtain for the absorption 
coefficient 



A 



1 



E- 



JE* 



E*^^\ 
dz J 



(70) 



likL \ Eq\^\ dz 

where ki^ is the wave number of the field in the surrounding medium to the left of the 
sample, and -Eq is the electric field amplitude of the incoming radiation. Multiplying Eq. (I7j) 
without sources and its conjugate by E* and E^ respectively, and integrating over the entire 
structure one obtains 



A 



E 



dz ^^{z)E{z] 



QW 



Im 



(71) 



This expression has clear physical meaning. It shows that the absorption of the resonant 
photonic crystal is the sum of independent contributions of all quantum wells. Each contri- 
bution has an expected form of the product of the imaginary part of the exciton susceptibility 
and a term proportional to the projection of the em field onto the exciton state in the well. 
In order to find a relation between the absorption and the emission spectra we calculate the 
absorption of the wave of unit intensity incident normally from the left. Electric field inside 
each well, and respective Integra 
technique (see, for instance Ref. 



s in Eq. fjTT]) . can be found using standard transfer-matrix 
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n 



As a result Eq. (I7T1) can be presented in the form 
explicitly containing the same Green's functions as one used in this paper, Eq. fHTl) : 

2 



A 



E 



Im Xr, 



9 Xr. 



(72) 



Comparing this expression with Eq. (HOl) (with terms proportional to Fi and F2 omitted) 
we can relate the emitted intensity to the absorption coefficient for a single m-th well of the 
structure 



X_('m) 



-A{m) 



A(^m) 



(73) 



c Imxm c 7 

where, in the second expression, we assumed that exciton susceptibility has a Lorentzian form 

and is described by Eq. (fT8l) . This assumption can be violated under several circumstances. 
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for instance, when several exciton levels spectrally overlap and contribute to the emission,— 
or in the presence of the inhomogeneous broadening of excitons. If the latter case is treated 
in the effective medium approximation, as was discussed previously in this paper and in 



Refs. 



3G 
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it results in an effective susceptibility with non-Lorentzian shape. 
In order to derive a global relation between emission and absorption for the entire struc- 
ture, Eq. (173!) has to be summed over all wells. If all wells are identical, i. e. the source 
functions and non-radiative decay rates in all wells are the same, the summation is trivial 
and the global relation between absorption and emission coefficients has again the form of 
Eq. ( 1731) . Comparing this result with Eq.( l68i) we can establish relationship between mi- 
croscopic parameters such as the strength of the exciton-light coupling, characterized by 
a, non-radiative decay rate 7, the source function H, which is proportional to polarization 
correlation function, and photon distribution function, A'^^^'': 



^Kki^) = - exp {-n/ksT)^^ (74) 



Since the photon distribution function is practically independent of frequency on the scale 
of frequencies considered in this paper, Eq. (17^ is consistent with an assumption that both 
7 and are frequency independent quantities. Taking into account also that l/u term in 
Eq. (l73l) also changes very weakly on the same scale, we can conclude on the basis of both 
Eq. (!67I) and Eq. (!73|) that emission and absorption spectra in our case are directly propor- 
tional to each other: X_(cj) oc A{u). Numerical evaluation of the respective expressions 
completely confirm this conclusion. 

The assumption of all wells being the same, however, may violate in photoluminescent 
experiments due to attenuation of the pumping radiation. As a result, different wells may 
be characterized by different exciton distribution functions, and consequently by different 
source functions. In this case, while Eq. (l73|) remains valid locally for any particular well, its 
global version is not true anymore. This results in loss of the direct proportionality between 
absorption and emission spectra of our structures. This is clearly seen from FigJTl which 
shows modification of emission intensity caused by attenuation of pump, while absorption 
spectra obviously are not affected by this circumstance. 
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V. CONCLUSION 



In the present paper we studied spectrum of non-coherent radiation emitted by one- 
dimensional resonant photonic crystal structures. While for concreteness we focused on ex- 
citon luminescence in multiple-quantum-well structures, the general theoretical framework 
developed in this work can be applied to other structures of this sort. The results obtained 
in the paper can be classified in two groups. First, we have developed a powerful method 
of solving general linear response type of problems for one-dimensional layered structures of 
general type. The problem of luminescence of one-dimensional resonant photonic structures 
is just one example of such problems, in which one is looking for the radiative response 
of the system caused by incoherent periodically distributed emitters. Our approach allows 
expressing Greens' function of the structure in terms of transfer matrices describing propa- 
gation of the radiation through the system. As a result we are able to present the spectrum 
of the emitted light in terms of reflection and transmission coefficients of the structure in 
question. Also, with the help of a special version of transfer-matrix description of transmis- 
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and 
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sion/refiection properties of resonant photonic crystals developed in recent Refs 
we were able to obtain a closed analytical expression for the spectrum of luminescence of 
an resonance photonic crystal structure with an arbitrary number of identical periods. An 
important characteristics of these general results is that they are obtained in terms of par- 
ticular solutions of an initial value problem for a structure with an arbitrary spatial profile 
of the refractive index. The latter problem can always be easily solved either analytically, 
or in most cases at least numerically, and, therefore, emission characteristics are expressed 
in our approach in terms of easily accessible quantities. 

The second group of the results is concerned with application of our general formalism to 
the particular case of Bragg or near-Bragg multiple-quantum-well structures. We analyzed 
the luminescence spectrum of these structures and established its main qualitative and 
quantitative characteristics. In particular, we explained the absence of luminescence in the 
spectral region of polariton stop-band, which was shown to be due to a combination of 
two factors: diminishing of transmission of light through the structure in the vicinity of the 
exciton frequency, and a significant spectral separation of the latter from transparent regions 
because of formation of the wide band-gap. It is interesting to note that this result agrees 



with quantum-field calculations of Ref. 



25, where similar effect was called " subradiance" . 
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Our calculations show, however, that this effect can be explained on purely classical ground. 

We also considered modification of the spectrum when the period of the structure becomes 
slightly de-tuned from the exact Bragg conditions. Comparison of our calculations with ex- 
perimental spectra demonstrated an important role played by inhomogeneous broadening of 
excitons in formation of the spectra of luminescence of the structures under consideration. 
We also showed that these spectra are influenced by a great deal of other effects such as 
attenuation of the pump, or presence of cladding layers in the structure. We found that our 
calculations produce good agreement with experiment in terms of positions of the peaks of 
the luminescence, but not for the relative height of the peaks. This, however, is not very 
surprising, because the intensities of the peaks depend on many various circumstances. One 
of the most intriguing effect, which was not considered in the paper, but which could signifi- 
cantly influence the distribution of the luminescence intensity between its peaks, is acoustic 
phonon-induced scattering between collective exciton polariton states formed in quasi-Bragg 
MQW structures. This possibility is supported by the fact that spectral separation between 
luminescence peaks in typical experiments^^"^ is of the same order of magnitude [ImeV) as 
an average energy of acoustic phonons in GaAs. Consideration of this effect is out of the 
scope of the current paper, but will be presented in the subsequent publications. 

In order to achieve a better understanding of luminescence of the quasi-Bragg structures, 
one needs additional experimental data, which would provide information for assessing the 
role of different effects. For instance, since one can, to some extent, control inhomogeneous 
broadening of excitons by growth conditions, luminescence measurements on a series on 
sample grown under different conditions could clarify the role of inhomogeneous broadening. 
Another possibility can be to excite luminescence by pumping the sample from both sides, 
which will reduce effects due to inhomogeneity of source function, and assess the role of this 
effect. Temperature dependence of the intensity of the luminescent maxima can be used to 
verify the role of acoustic phonon scattering. 

An interested question studied in the last section of the paper is concerned with re- 
lation between luminescence and absorption spectra in multiple-quantum-well structures. 
First, using thermodynamical arguments we derived a version of Kirchhoff's law specifi- 
cally adapted for one-dimensional structures under consideration. Then, using developed 
formalism, we were able to establish the relation between luminescence and absorption in 
terms of microscopic characteristics of the system such as polarization correlation function 
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and exciton susceptibility. Comparing the two results we established a relation between the 
microscopic parameters consistent with the Kirchhoff's law. In particular, we found that 
if the spectral region of interest is small compared to the characteristic energy scale of the 
photon distribution function, both polarization correlation function and the non-radiative 
decay rate of excitons can be considered as frequency independent. 
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Appendix 

In this Appendix we will provide a sketch of quantum-mechanical calculations of lumines- 
cence from a single quantum well. The objective of this exercise is to provide a microscopic 
justification for the quasi-classical approach employed in the paper, and demonstrate rela- 
tions between phenomeno logical source function S(p, t) and microscopical characteristics of 
the system. This relation can be derived only by using the second quantization of the exciton 
and photon states and a relevant microscopic kinetic equation. Here we will demonstrate 
the second-quantization approach for a single QW-structure with a QW layer sandwiched 
between semi-infinite barriers. For simplicity, the dielectric contrast — is set to zero. 

The 2D -space is naturally divided into two regions, radiative and nonradiative, respec- 
tively, with k < ko = {uJo/c)nb and k > k^. In the process of relaxation of electron-hole 
pairs non-radiative excitons with k > ko are being populated first^, and we assume that, 
for these excitons the criterion kl ^ 1 for validity of the Boltzmann kinetic equation is 
fulfilled. Here k'^ is the average value of k'^ describing the spread of exciton population in 
the fc-space, / is the 2D exciton mean free-path length, / = {hk/M)Tp, Tp is the exciton 
momentum scattering time, and M is the exciton in-plane effective mass. For 3D photons 
we introduce a quantization box of the volume V = SL with the macroscopic interface area 
S and the macroscopic length L along the QW growth direction. 
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In order to calculate the spectral intensity of light emitted by excitons in the QW we 
apply the Keldysh diagram technique^i^, in a way similar to one used for the description of 
exciton or exciton-polariton photoluminescence in bulk crystals, see for instance, Ref. |55|]56 . 
In this method the intensity can be written as 

Ig = hUgWg , Wg= - ^Ihll^ j D^^ (1^ . (75) 

Here Wg is the emission rate of a photon with the 3D photon wave vector g, is the 

photon Green function presented in the left-hand side of the diagram equation shown in 
Fig. [11] by an short-dashed line connecting the lower and higher horizontal parts of the 
Keldysh contour labeled "+" and respectively; 7 is a positive photon damping rate 
which is introduced for the formal reasons in order to stabilize the photon distribution in 
the steady-state regime [(27) is the photon lifetime in the quantization box]. In the end 
7 is set to +0 since in a single QW structure with the microscopic length L the actual 
photon damping rate is negligible. Equation (1751) is derived taking into account that, in the 
steady-state regime, one has 

1 7 

-n {uj - uOgY + Y 



Dg: = -2nNg „ ) , (76) 



where Uq = cq/rif,, Ng = Wg/{2'y) is the steady-state photon distribution function. Here 
the temperature is assumed to be very small as compared to the exciton excitation energy, 
which allows one to neglect the equilibrium photon-state population. Note that the close-to- 
normal energy flux in the frequency range duj and within the area (Pq\\ in the plane [qx, Qy) 
is related with Ig by 

d^q\\duj Lnb 

The Green function D~'^ is found from the diagram equation presented in Fig. [TTJ The 
diagrams on the right-hand side of the equation describe the acoustic-phonon assisted scat- 
tering of an exciton from the non-radiant state k' to the radiative state gy, where q\\ is the 
in-plane component of g; Q and uoq are the phonon wave vector and frequency. In Fig. [11] 
the short-dashed lines mean i-D*^, where the superscripts s', s = ± show the position of the 
ends of the photon Green functions, 

D— = - = . 
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FIG. 11: . 

The solid lines represent exciton Green functions iCfj, given by 

1 



n — ^ _ ic^+V = 

^(exc)^-p^ and Fofc are the exciton excitation energy, non-radiative and radiative damping 
rates, fk' is the exciton distribution function assumed to be small, fk> <^ 1- Each vortex 
represents the matrix element of exciton-phonon (Vg) or exciton-photon {Mq^^) interaction 
multiplied by ±i/^, plus for the lower part and minus for the upper part of the Keldysh 
contour. Since we neglect the damping of acoustic phonons the upward and downward 
phonon (long-dashed) lines can simply be replaced by the factors tuq and rriQ + 1, where rriQ 
is the phonon occupation number. For a process described by any vortex the conservation 
of the in-plane wave vectors should be satisfied. In particular, k' = k ± Q\\ for the phonon 
emission and absorption, respectively. As a result, we obtain 



k',Q 

where summation over k/ is carried only over non-radiative states because the population of 
the radiative states is assumed to be negligible. It follows then that the photon generation 
rate Wg in Eq. (175!) can be presented in the form 

27r, _ ,2 1 Ffc + Fofc 



^{iu-Jr^y + iTk + Tokr 
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k',Q 

+ {rnq + 1) 6{uj + - c^i^'"'') 4',fc+Q||] ■ 

Here, taking into account the condition kl ^ 1 we replaced in the last equation the phonon's 
resonance Lorentzians by Dirac delta-functions. 

In the Bragg and quasi-Bragg structures, the radiative excitonic states \n, k) in different 
wells n with the same wave vector k are strongly coupled by the electromagnetic field. In 
contrast, excitons with k satisfying the conditions k > k^, a/A;^ — k'^d ^ 1 are decoupled and 
can be considered as excitations isolated in particular wells. The photon emission outgoing 
from a given QW can be described by Eq. ( 1771) as well or, equivalently. On the other hand, 
the propagation of this outgoing light wave through the whole MQW structure accompanied 
by its reabsorption and escape into the vacuum or substrate, can be described classically. 

Instead of the above microscopical consideration, in the main part of the paper the 
photoluminescence spectra are calculated with the help of a quasi-classical Langevin-like 
approach by introducing the random source uj) into expression for exciton polarization, 

see Eq. (jS]). Now, the explicit expression for the correlator Smiw, k) of the random sources 
in Eq. ( H6|) can be readily found from Eq. ( 1771) . In particular, one has 

(78) 

where Tq is the radiative damping rate of an exciton with k = excited in a single QW 
structure. 

The applicability of the Langevin method can be justified by considering Heisenberg 
equations of motion for exciton annihilation (cfc) and creation (c^) operators. For the exciton- 
phonon Hamiltonian 

kk'Q 

where dgydq are the acoustic-phonon annihilation and creation operators, one obtains 

dk{t) = -Jr^Ckit) - ^Mfbkit) + ikit) (79) 
with bk being the photon annihilation operator. Here term 

ikit) = ^ (VQdQ6k',k~Q^^ + VQdQ^k',k+Q^^ Ck' 

k'Q 
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2nq 



^m{z)dz 



playes the role of the Langevin source operator, and summation again is carried over only 
values of kl corresponding to non-radiative states. Fourier transform of the correlator of this 
operator 



is equal to W^^^°'^\uj,k), which agrees with Eq. (ITS]) if we take into account the relation 
between the 2D exciton envelope function and exciton-induced polarization Pexc- 
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